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Abstract 

We extend the complete Mellin (CM) representation of Feynman am- 
plitudes to the non-commutative quantum field theories. This represen- 
tation is a versatile tool. It provides a quick proof of meromorphy of 
Feynman amplitudes in parameters such as the dimension of space-time. 
In particular it paves the road for the dimensional renormalization of 
these theories. This complete Mellin representation also allows the study 
of asymptotic behavior under rescaling of arbitrary subsets of external 
invariants of any Feynman amplitude. 

1 Introduction 

Recently non-commutative quantum field theories such as the and the Gross- 
Neveu2 models have been shown renormalizable [Ij El [31 El [5l |6] provided the 
propagator is modified to obey Langmann-Szabo duality [7] . These theories are 
called "vulcanized" , and in [51 IH] the Feynman-Schwinger parametric represen- 
tation was extended to them. For recent reviews, see [TUl lllj . 

In this paper we perform a further step, generalizing to these vulcanized 
NCQFT the complete Mellin (CM) representation. This CM representation is 
derived from the parametric representation established in [8l [9] and provides a 
starting point for the study of dimensional renormalization and of the asymp- 
totic behavior of Feynman amplitudes under arbitrary rescaling of external in- 
variants. 

The study of asymptotic behaviors in conventional (commutative) field theo- 
ries started in the 1970's with the papers [12l[T3l[T4] which use the BPHZ renor- 



malization scheme. Their approach is based on the Feynman-Schwinger para- 
metric representation of amplitudes. This representation involves "Symanzik" 
or "topological" polynomials associated to the Feynman diagram. In commu- 
tative theories these "topological polynomials" are written as sums over the 
spanning trees or 2-trees of the diagram [Tallin]. The corresponding mathemat- 
ical theory goes back to the famous tree matrix theorem of Birkhoff (see [T7] 
and references therein). 

In vulcanized NCQFTs, propagators are no longer based on the heat kernel 
but on the Mehler kernel. The kernel being still quadratic in position space, 
explicit integration over all space variables is still possible. It leads to the NC- 
QFT parametric representation. This representation no longer involves ordinary 
polynomials in the Schwinger parameters, but hyperbolic polynomials [S]. As 
diagrams in NCQFT are ribbon graphs, these hyperbolic polynomials contain 
richer topological information than in the commutative case. In particular they 
depend on the genus of the Riemann surface on which the graphs are defined. 

Returning to the Mellin transform technique, it was introduced in commu- 
tative field theory to prove theorems on the asymptotic expansion of Feynman 
diagrams. The general idea was to prove the existence of an asymptotic series 
in powers of A and powers of logarithms of A, under rescaling by A of some 
(Euclidean) external invariants associated to the Feynman amplitudes. The 
Mellin transform with respect to the scaling parameter allows to obtain this re- 
sult in some cases and to compute the series coefficients, even for renormalized 
diagrams |18j . But it did not work for arbitrary subsets of invariants. 

In the subsequent years, the subject advanced further. The rescaling of inter- 
nal squared masses (in order to study the infrared behavior of the amplitudes) 
was treated in |^. There the concept of "FINE" polynomials was introduced 
(that is, those being factorizable in each Hepp sector ^\ of the variable^. 
It was then argued that the Mellin transform may be "desingularized" , which 
means that the integrand of the inverse Mellin transform (which gives back the 
Feynman amplitude as a function of A) has a meromorphic structure, so that the 
residues of its various poles generate the asymptotic expansion in A. However, 
this is not the case under arbitrary rescaling, because in many diagrams the 
FINE property simply does not occur. 

A first solution to this problem was presented in jSUj by introducing the so- 
called "multiple Mellin" representation, which consists in splitting the Symanzik 
polynomials in a certain number of pieces, each one of which having the FINE 
property. Then, after scaling by the parameter A, an asymptotic expansion can 
be obtained as a sum over all Hepp sectors. This is always possible if one adopts, 
as done in [2TJ [231 US] , the extreme point of view to split the Symanzik polyno- 
mials in all its monomials. Moreover, this apparent complication is compensated 
by the fact that one can dispense with the use of Hepp sectors altogether. This 
is called the "complete Mellin" (CM) representation. 

The CM representation provides a general proof of the existence of an asymp- 
totic expansion in powers and powers of logarithms of the scaling parameter in 

Hepp sector is a complete ordering of the Schwinger parameters. 
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the most general case. Moreover the integrations over the Schwinger parameters 
can be exphcitly performed, and we are left with the pure geometrical study of 
convex polyhedra in the Mellin variables. The results of [20] are obtained in a 
simpler way |21j . and asymptotic expansions are computed in a more compact 
form, without any division of the integral into Hepp sectors. 

Moreover the CM representation allows a unified treatment of the asymptotic 
behavior of both ultraviolet convergent and divergent renormalized amplitudes. 
Indeed, as shown in [SUES], the renormalization procedure does not alter the 
algebraic structure of integrands in the CM representation. It only changes the 
set of relevant integration domains in the Mellin variables. The method allows 
the study of dimensional regularization [211 US] and of the infrared behavior of 
amplitudes relevant to critical phenomena |24| . 

The CM representation is up to now the only tool which provides asymptotic 
expansions of Feynman amplitudes in powers and powers of logarithms in the 
most general rescaling regimes. More precisely consider a Feynman amplitude 
G(sfc) written in terms of its external invariants Sfc (including particle masses), 
and an asymptotic regime defined by 

Sfe X'"'Sk, (1) 

where at may be positive, negative, or zero. Let A go to infinity (in this way 
both ultraviolet and infrared behaviors are treated on the same footing) . Then 
G{sk) as a function of A has an asymptotic expansion of the form: 

GiX;sk)= E GMXnii'^X, (2) 

P=Pmax 9=0 

where p runs over decreasing rational values, with Pmax as leading power, and 
q, for a given p, runs over a finite set of non-negative integer values. 

The CM representation proves this result for any commutative field theory, 
including gauge theories in arbitrary gauges [25^ . 

The above general result should hold mutatis mutandis for non-commutative 
field theories. In this paper we provide a starting point for this analysis by 



constructing the CM representation for cj) (Theorem 4.1 below). Moreover 



using this representation we check meromorphy of the Feynman amplitudes in 



the dimension of space-time (Theorem 5.1 1 



The main difference with the commutative case is that this integral repre- 
sentation (previously true in the sense of functions of the external invariants) 
now holds only in the sense of distributions. Indeed the distributional character 
of commutative amplitudes (in momentum space) reduces to a single overall 
5-function of momentum conservation. This is no longer true for vulcanized 
NCQFT amplitudes, which must be interpreted as distributions to be smeared 
against test functions of the external variables. 
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2 Commutative CM representation 

For simplicity we consider a scalar Feynman amplitude. To get the CM repre- 
sentation [21\ . rewrite the Symanzik polynomials as 



= E n ^ E ' ^(") = E f n ^ E 

j i=\ j k \£=1 / k 



(3) 



J «=1 3 

where j runs over the set of spanning trees and k over the set of the 2-trees, 

_ J if the line £ belongs to the 1-tree j , . 

"^■^ ~ [ 1 otherwise ^ ' 

and 

r if the line I belongs to the 2-tree k 

''^'^ = I 1 otherwise. 

The Mellin transform relies on the fact that for any function f{u), piecewise 
smooth for u > 0, if the integral 

g{x)= / duu-'^-'fiu) (6) 
Jo 

is absolutely convergent for a < Re x < (3, then for a < a < (3 



1 



(T+'ioo 



f{u)^T^. dxg{x)u\ (7) 

Consider D the space-time dimension to be for the moment real positive. 
Taking f{u) = e~", and applying to m = Vk/U one gets 

e"^^/^ = / T{-y,){^y\ (8) 

where J^^ is a short notation for d{im^k) ^ ^--^^Yi Re fixed at r^- < 0. We 
may now recall the identity 

r(u)(^ + B)""= / "^P"^ ''\ {~x)A'V{x + (9) 
J-co 27r 

Taking A = ?7i(x) and B = U2 + U3 + ■ ■ ■ and using iteratively the identity 
above, leads, for u = yk + D/2, to 



r(E2/^+?) c/-^^--*= /nr(-.)c^r. 



(10) 



with Re Xj — aj < 0, Re (X^fc 2/fc + y) = X^fe''"*; + T > 0' ™d / means 
I-^Uj ^^^^ with J:^Xj + Efeyfe = -f • Then, using (g and the 
amplitude is written as 



TT Y(—x ) r°° 
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where 

(j)i = ^UijXj +^Vikyk + l- (12) 

j k 

The symbol means integration over the independent variables ^"^^^ , 

in the convex domain A defined by (cr and r standing respectively for Re Xj 

and Re yk): 



A = ->a,T 



CTj < 0; Tfc < 0; J2j ^] + Efc Vk = ; 
Vi, Re (f>i = J2j UijCTj + J2k ^ikTk + 1 > 



(13) 



This domain A is non-empty as long as D is positive and small enough so that 
every subgraph of G has convergent power counting |21j , hence in particular for 
the (j)* theory it is always non empty for any graph for < I? < 2. 

The a integrations may be performed, using the well-known representation 
for the gamma function, so that we have 



ii e 



^'ap-'^T{c^OHr'" (14) 



and we finally get the CM representation of the amplitude in the scalar case: 
Ia{sk,mf)^ [ ^^^^^T[sfT{^y,)l[{rnf)-^^T{^,). (15) 

This representation can now be extended to complex values of D. For instance 
for a massive cj)'^ graph, it is analytic in D for < 2, and meromorphic in D 
in the whole complex plane with singularities at rational values; furthermore its 
dimensional analytic continuation has the same unchanged CM integrand but 
translated integration contours ^22, |23j . 



3 Non- commutative parametric representation 

Let us summarize the results of Define the antisymmetric matrix a as 

The (5— functions appearing in the vertices can be rewritten as an integral over 
some new variables pv called hypermomenta, via the relation 

(17) 

Let 5 be a ribbon graph. Choosing a particular root vertex V , the parametric 
representation for the amplitude of Q is expressed in terms of tg = tanhaf/2 
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(where ae is the former Schwinger parameter) as 

Ag^ J^'lldt,{l-tj)^-' J da;dpe-9^'=^* (18) 

where D is the space-time dimension, X summarizes all positions and hyper- 
momenta and G is a certain quadratic form. CaUing and py the external 
variables and Xi,pi the internal ones, we decompose G into an internal quadratic 
form Q, an external one M and a coupling part P so that 

X = {xe Pv Xi p^) , G=lpt , (19) 

Performing the Gaussian integration over all internal variables one gets the 
non-commutative parametric representation: 

This representation leads to new polynomials HUg y and HVg y in the ti 
{£ = 1, . . . ,L) variables, analogs of the Symanzik polynomials U and V of the 
commutative case, through 



Ag{{xe}, Pv) ^K' f Y[dU{l-tl)'i- 



HUgy{t) 



(21) 



The main results of O |9] are 



• The polynomial HUg y and the real part of HVg y have a positivity prop- 
erty. They are sums of monomials with positive integer coefficients, which 
are squares of Pfaffians with integer entries. 

• Leading terms can be identified in a given "Hepp sector" , at least for 
orientable graphs. In HUg y they correspond to hyper-trees which are the 
disjoint union of a tree in the direct graph and an other tree in the dual 
graph. Any connected graph has such hypertrees. Similarly "hyper-two- 
trees" govern the leading behavior of HVg y in any Hepp sector. 

Let us relabel the 2L internal positions of as L short and L long variables in 
the terminology of [5J . It has been shown in [8 that for the Grosse-Wulkenhaar 
(j)*'^ model: 

HUg.yit) = S'""'"- '^LH*^ n^^' 

Ku=IUJ, n+\Ku \ odd 1<^I i'£J 

= E«A-.nc'^'^E^t^^- (22) 

Ku I Ku 

where 
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• / is a subset of the first L indices (corresponding to the short variables) 
with |/| elements, and J a subset of the next L indices (corresponding to 
the long ones) with | J| elements. 

• S is the antisymmetric part of the quadratic form Q restricted to these 
2L short and long variables (i.e. omitting hypermomenta) 

• riKu = Pf(-B^) is the Pfafhan of the antisymmetric matrix obtained from 
B by deleting the lines and columns in the set Ku = lU J. 

• kx^, is \Kjj\ — L — F + 1, where F is the number of faces of the graph, g 

is the genus of Q and s is a constant. 

• = s^^~'"'un%^, and 

To if £ e / and £ J 
UiKu = { 1 if (£ ^ / and £ ^ J) or (£ e / and £ ^ J) . (23) 
[2 if ^ ^ 7 and ^ e J 

The second polynomial HV has both a real part HV^ and an imaginary 
part HV^ , more difficult to write down. We need to introduce beyond / and J 
as above a particular line r ^ / which is the analog of a two-tree cut. We define 
Pf (i?^^ .) as the Pfaffian of the matrix obtained from B by deleting the lines and 
columns in the sets /, J and r. Moreover we define e/^^ to be the signature of 
the permutation obtained from (1, . . . , c?) by extracting the positions belonging 
to / and replacing them at the end in the order 



• • • • • • )«- 



• • ■,d,ir,i\i\ . . . ,«! . 



(24) 



where d is the dimension of the matrix Q. Then: 





(25) 



where 




2 



E^^ E 



er' 



£K.xPf(S^Vf) 



(26) 



and ViKv given by the same formula as u^k^ 
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The imaginary part involves pairs of lines t, t' and corresponding signatures 
(see dlin] for details): 



Kv=iuJe^i £'eJ 



ei,e2 rr' 



where 

4. - eK.Pf(i?^v) f E(E ^er£K.rr'Pf(SKVfr')^e'r')a;e^a:,, ) . (28) 



e.e' r,r' 



A similar development exists for the LSZ model (^), although more in- 
volved. The first polynomial in this case writes 

ic{i...L}, lei ^ i'e{i,...,L} 

n+\I\ odd 



E aK^Wer- ^Y.HUku- (29) 



Ku=IU.J I Ku 

Again / runs over the subsets of {1, ... , L}. But J now runs over subsets of / 
(representing the expansion of H^e/Il + ^^l); has again 

To if £ e / and £ e J , 
utKu^l I mil , aKu = T^^s'^a-k'n'i. (30) 

[2 if £ e / and £ ^ J ^ 

There are similar expressions for the second polynomial. 

To summarize, the main differences of the NC parametric representation 
with respect to the commutative case are: 



• the presence of the constants Oj in the form (22 1 of HU , 



• the presence of an imaginary part iHV in HV , see (27) 



• the fact that the parameters u^j and vg^ in the formulas above can have 
also the value 2 (and not only or 1). 
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4 Non- commutative CM representation 

For the real part HV^ of HV one uses again the identity (jsl as 



BY, 



U 



(31) 



which introduces the set of MeUin parameters . 

However for the imaginary part one cannot apply anymore the same identity. 
Nevertheless it remains true in the sense of distributions. More precisely 
we have for HV^^/U > and -1 < r^^ < 



n-VKv) 



U 



(32) 



The proof is given in Appendix B. 

Note that this introduces another set of Mellin parameters . The dis- 
tributional sense of formula ( 32 ) is a translation of the distributional form of a 



vertex contribution, and the major difference with respect to the commutative 
case. 



For the polynomial HU one can use again the formula ( 10 1 which rewrites 
here as 



\Kv 



D 



(HU) 



* = / (33) 



Note that the role of the Mellin parameters i/Kv of the commutative case is now 
played by the sum y'^^ + yj^^ . 



As in the commutative case, we now insert the distribution formulas (|31l 
; general form of the Feynman am 



(32) and (33 1 in the general form of the Feynman amplitude. This gives 



An = K 



where 



(34) 



Kv 



Here means integration over the variables — 



Im XK.- Im 



JA 

A is the convex domain: 



i ' 27ri 



and 



^Ku < 0; r^,, < 0; -1< rj^ < 0; 

Y.Ku + UkAVkv +yKv)^- 
Re (jji = J2ku '^^Ku^Ku 

+ EKA^^Kvr^v + ^iKv^D + 1 > 



D. 

2 ' 



(35) 

2^, , where 
(36) 
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and cr, t'^ and r-^ stand for Re xku^ Re and Re j/J 



The dti integrations in (34 1 may be performed using the representation for 
the beta function 

ndi,(l-t?)f-4^^ = l/3(^|,|). (37) 



Furthermore one has 



r(f)r(f) 



2 ' 2 r(^^±^) 

This representation is convergent for < < 2 and we get: 

Theorem 4.1. Any Feynman amplitude of a ip** graph is analytic at least in 
the strip < < 2 where it admits the following CM representation 



which holds as tempered distribution of the external invariants. 



(38) 



Proof We have to check that for any (f>** graph, the domain A is non empty for 
< < 2. It is obvious to check that A is non empty for < JRD < 1; indeed 
since the integers u and v are bounded by 2, we can put an arbitrary single 
xku close to —D/2 and the others with negative real parts close to and the 
conditions (pi > will be satisfied. With a little extra care, one also gets easily 
that A is non empty for < < 2. Indeed one can put again one single xk^ 
close to —D/2, provided for any £ we have u^^Ku 2- For the ordinary Grosse- 
Wulkenhaar model one can always find Ku — I U J where / — {!,..., L} is the 



full set of lines 0, hence u^^Ku 2 for any ^ by (23 1. In the "covariant" case 
treated in [2], one can simply take Ku with J — I. Considering (30) ensures 
again obviously that ui^k^ 7^ 2 for any £. 

The analyticity statement is part of a more precise meromorphy statement 
proved in the next section. □ 

We have thus obtained the complete Mellin representation of Feynman 
amplitudes for non-commutative QFT. Let us point out first that the 
distributional aspect does not affect the last product, the one involving the 
dimensionality of space-time D. Therefore, one can still prove the meromorphy 
of the Feynman amplitude Ag{D) in the space-time dimension D essentially as 
in the commutative case. This is done in the next section. 



Another difference is that the dt( integration brought /3 functions (see (37l) 



hence more complicated ratios of F functions than in the commutative case. In 
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particular singularities appear for negative D which come from the integration 
near ti — 1, hence the "infrared side". They have absolutely no analog in the 
commutative case and are due to the hyperbolic nature of the Mehler kernel. 
At zero mass these singularities occur for D real negative even integer and in 
the massive case they occur for D real but sufficiently negative depending on 
the mass. 



5 Meromorphy in the space-time dimension D 



The complete Mellin representation (38 1 for non-commutative theories allows a 



quick proof of the meromorphy of the Feynman amplitude Ag in the variable 
D, since D appears only as argument of T functions which are meromorphic. 

Theorem 5.1. Any amplitude Ag is a tempered meromorphic distribution in D, 
that is Ag smeared against any fixed Schwarz-class test function of the external 
invariants yields a meromorphic function in D in the entire complex plane, with 
singularities located among a discrete rational set Sq which depends only on the 
graph G not on the test function. In the (f>* case, no such singularity can occur 
in the strip < < 2, a region which is therefore a germ of analyticity 

common to the whole theory. 

Proof r functions and their inverse are meromorphic in the entire complex 



plane. Formula (38 1 is a tempered distribution which is the Fourier transform of 
a slow-growth function, see Appendix A. It must be smeared against a Schwarz- 
class test function and by definition the result can be computed in Fourier space, 
as an ordinary integral of the product of the slow-growth function by the rapid 
decaying Fourier transform of the test function. At this stage the /^-dependence 
occurs in fact in the integrand of an ordinary integral. 

We can now copy the analysis of ^22l [23] . We introduce general notations 
Zi, for the variables x and y and the lattices of polar varieties — — where 
the different F functions have their singularities. This defines convex cells, 
and singularities in D can occur only when the hyperplane P{D) = {X^jy = 
—D/2} crosses the vertices of that lattice. Because the coefficients u and v can 
only take integer values (in fact 0, 1 and 2) this give a discrete set of rational 
values Sq which depend on the graph, but not on the test-function. Out of these 
values, to check that this integral is holomorphic around any dimension which 
does not cross the lattice vertices is a simple consequence of commuting the 
Cauchy-Riemann djy operator with the integral through Lebesgue's dominated 
convergence theorem. 

Finally to check meromorphy, fix a particular Dq S Sq. We can write 
the Laurent series expansion of the integrand near that ZJq which starts as 
ano{D - fo)"""- Multiplying by (Z) - Dq)"" we get a finite limit as D ^ Do, 
hence the isolated singularity at Do cannot be essential. Repeating the argument 
at each point of Sg we obtain meromorphy in the space time dimension D in 
the whole complex plane C, apart from some discrete rational set where poles 
of finite order can occur. □ 
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Recall that this theorem is the starting point for the dimensional renormal- 
ization of NCQFTs which will be published elsewhere [25] . 



Appendix A: Proof of the formula (32 1 



Let us return to formula ( 32 ) . It follows from 



+00 7, 

X{u)^x{u)l r(-,s-zt)(zur+'V- (39) 



Lemma 5.1. For — 1 < s < 0, as tempered distribution in u 

dt 
2n 

where x is the Heaviside step function ensuring u > 0. 

Proof The left hand side is a locally integrable function bounded by one (to- 
gether with all its derivatives) times the Heaviside function, so it is a tempered 
distribution. The right hand side is in fact a Fourier transform in t but taken at 
a value logu. To check that this is the same tempered distribution as e~'"x('^'); 
let us apply it to a Schwartz test function /. Writing the usual representation 
for the r function we have to compute 

I{f) = — / f{u)du / {iuy+'*dt / x-'-'^-^e-'^dx. (40) 

271' Jo J-oo Jo 

But the integrand is analytic in x in the open upper right quarter of the complex 
plane. Let us rotate the contour from x to ix. The contribution of the quarter 
circle with radius sent to infinity tends to zero because of the condition — 1 < 
s < 0, and we find 

-t-oo 



1 poo r-\-QO Pit 

I{f) = — / f{u)du / {iuY+^'dt lim / [ix)-' 
■^T^ Jo J-oo R^°°Jo 

1 r+oo /•oo /•/?, 

= dt f{e'")e"^''+^^e'^'"dv lim / x-'-'^-^e-'^'dx 

/+00 /^R 
dtG{t) lim / x-'-'^-^e'^'^'dx 
-oo Jo 



1 

2tt 



pR 1 p+oo 

lim / dxe-'^'x-"-^— / dtG(i)e-*"°s^ 
Jo 27r J_^ 

f{x)e-~"dx . (41) 

where G{v) — f[e^)e^'^^^^'> is still rapid decay and in the second line we change 
variables from u to v — log u and use the definition of Fourier transform for 
distributions. □ 

Appendix B: The massive case 



In the massive case we have to modify equation ( 37 ) to include the additional 
-aem pgj, unfortunately we have e~"* = so that (37) becomes 

/'n*(i-'?)*-'(^)'""'f'-=ift4'?'- « 
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where the modified "/Sm" function is no longer an explicit quotient of F functions, 
and depends of the mass m. 

It is rather easy to check that the analyticity and meromorphy properties of 
this modified beta function are similar to those of the ordinary beta function, 
at least for "SiD > 0. Indeed 

ft4f.f)-«f,f)^2/n.,(l-.?)*-'[(i^)"""-ll<?-M43) 

The integral above is convergent for tg close to 1 if 'StD > and for ti close 
to if 3?(/); > — 1 so that it is an analytic function of D and (f>i in this domain, 
hence /3 and /3m have the same singularity structure. 

Nevertheless this CM representation becomes less explicit and therefore less 
attractive in this massive case. Remark however that for > 0, masses 
are not essential to the analysis of vulcanized NC field theories which have no 
infrared divergencies and only "half-a-direction" for their renormalization group 
anyway. As remarked earlier, masses can only push further away the infrared 
singularities at negative D that come form the hyperbolic nature of the Mehler 
kernel. But it is not clear whether such infrared singularities at negative non 
commutative dimensions have any physical interpretation. 
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